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ABSTRACT 

The study of time lags between spikes in Gamma-Ray Burst light curves in 
different energy bands as a function of redshift may lead to the detection of effects 
due to Quantum Gravity. We present an analysis of 15 Gamma-Ray Bursts with 
measured redshift, detected by the HETE-2 mission in order to measure time 
lags related to astrophysical effects and search for Quantum Gravity signature in 
the framework of an extra-dimension string model. 

The wavelet transform method is used both for de-noising the light curves 
and for the detection of sharp transitions. The use of photon tagged data allows 
us to consider various energy ranges and to evaluate systematic effects due to 
selections and cuts. 

The analysis of maxima and minima of the light curves leads to no significant 
Quantum Gravity effect. A lower limit at 95% Confidence Level on the Quantum 
Gravity scale parameter of 2 x 10^^ GeV is set. 

Subject headings: gamma rays: bursts — time lags — wavelet analysis — quan- 
tum gravity 



Introduction 



Particle Physics pr ovides a challenging des c ription of Quantum Gravity in the framework 
of the String Theory JeUIs et all IiQQtI . 1 19981 : lAmelino-Camelia et al. l ll997[ l. Within this 
scheme, gravitation is considered as a gauge interaction and Quantum Gravity effects result 
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from graviton-like exchange in a background classical space-time. This approach does not 
imply a "spontaneous" Lorentz symmetry breaking, as it may appear in General Relativity 
with Loop Quantum Gravity, which postulates discrete space-t ime in the Planckian regime 
dGambini fc PuUinlllQQQl : ISmolinlboOll : klfaro et aliboool . l2002h . 



In recent years, several experimental probes have been propos ed to test Lorent z in- 
variance in t he frame of both particle physics and astrophysics (see iMattinglyl (120051 ) and 
Sarkarl (120021) for a revi e w). I n the domain of gamma-ray astronomy, the idea proposed by 
Amelino-Camelia et al.l (119981 ) to use Gamma-Ray Bursts (GRBs) to measure arrival time 
of photo ns of dif f erent energies was taken further and applied to other sources like pulsars or 
blazars. iKaaretl (119991 ) gets a limit on the q uantum gravity en ergy scale of 1.8 x 10^^ GeV 
using the Crab pulsar observed by EGRET. iBiller et al.l (119991 ) use a gamma-ray flare from 
Markarian 421 and obtain a limit of 6 x 10^^ GeV. 

GRBs are the most distant variable sources detected by present experiments in the 
energy range from keV to GeV. These violent and explosive events are followed by a delayed 
emission (an afterglow) at radio, infrared, visible and X-ray wavelengths. The energy released 
during the explosion phase is of the order of 10^^ erg when the beaming corrections are 
applied. 

As mentioned above, GRBs have been proposed to study modification of photon prop- 
agation implied by quantum gravity. They are good candidates for this kind of work since 
they are bright transient sources located at cosmologi cal distances . Some studies use only 
one GRB to get a limit on the quantum energy scale. iBoggs et al.l (120041 ) obtain a limit of 
1.8 X 10^^ GeV using GRB 021206 observed by RHESSI. This is currently the best limit 
obtained with a GRB but it reli es on an e stima tion of the redshift which suffers from a big 
uncertainty. With GRB 9301 31. ISchaefeil (Il999h gets a limit of 8.3 x 10^^ GeV with 30 keV 
to 80 MeV photons. Finally, IeUIs et al.l J2003I . bood ) make use of several GRBs at differ- 
ent distances with measured redshifts. Using 9 GRBs observed by BATSE and OSSE and 
35 GRBs seen by BATSE, HE TE-2 and S WIFT , they obtain limits which are respectively 
6 X 10^^ GeV and 9 x 10^^ GeV. IeUIs et al.l tood ) use public data of HETE, which is available 
in 3 fixed energy bands. In the present paper, we use photon tagged data, which allow us 
to study several ene rgy gap scenario s. All papers quoted in this paragraph make use of the 
model described in (lEUis et al.ll2000l) . However, a different formalism has been proposed by 
Rodriguez Martmez et al.l (120061 ) which leads to a limit of 6.6 x 10^^ GeV, after analysis of 
GRB 051221A. 

While GRBs constitute interesting sources to test the Lorentz violation, they are not 
perfect signals. The GRB prompt emission extends over many decades in energy (from the 
optical to GeV) and it is conceivable that the emission at very different wavelengths (e.g. 
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optical and gamma-rays) is produced by different mechanisms, resulting in different light 
curves. Until we fully understand the radiative transfer in GRBs, we must restrict ourselves 
to an energy domain where the emission is produced by a single process. This is the case for 
the prompt X-ray and low energy gamma-ray emissions which have very similar light curves 
and are thus appropriate for the study of Lorentz violation. The influence of the source 
effects in the present analysis will be addressed later in this article. 

In order to search for Quantum Gravity using the GRBs detected by the HETE-2 
mission ^ we concentrate on the extra-dimension string model interpretation. The employed 
model flEUis et al.l [2OO0I ) relies on the assumption that photons propagate in the vacuum 



which may exhibit a non-trivial refractive index due to its foamy structure on a characteristic 
scale approaching the Planck length / ~ "^piinck- "^^^^ would imply light velocity variation 
as a function of the energy of the photon (E). In particular, the effects of Quantum Gravity 
on the light group velocity v would lead to: 



v{E) = -—, (1) 

where n{E) is the refractive index of the foam. Generally, the Quantum Gravity energy scale 
Eqg is considered to be close to the Planck scale. This allows us to represent the standard 
photon dispersion relation with E/Eqg expansion: 

.2 5 I ^ - E ^, E'^ A , ^. . E 



cy = E'\l + ^-— + Oi^)\ and v{E)^c(l-^^], (2) 



where ^ is a model parameter whose value is set to 1 in the following (lAmelino-Camelia et al. 



19981). 



The analysis of time lags as a function of redshift requires a correction due to the 
expansion of the Universe, which depends on the cosmological model. Following t he analysis 
of th e BAT SE data and more recently of the HETE-2 and S WIFT GRB data JeUIs et al. 
2003 . 2006 ) and considering the Standard Cosmological Model ( Bahcall et al. 19991 ) with flat 



expanding Universe and a cosmological constant, the difference in the arrival time of two 
photons with energy difference AE is given by the formula: 



where 



(4) 
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We assume iltot = + = 1, = 0.7 and Hq = 71 km s ^ Mpc 



Equation [3] is different from the one used in lEUis et al.l (l2003l . |2006| ) where the authors 
computed the time lag at redshift z. As the time lag is measured at redshift 0, a factor {1 + z) 
has been introduced in the integral. A more detailed computation leading to Equation [3] is 
given in Appendix |Al 

In order to probe the energy dependence of the velocity of light induced by Quantum 
Gravity, we analyse the time lags as a function of the redshift. Possible effects intrinsic to 
the astrophysical sources could also produce energy dependent time lags. The analysis as a 
function of z ensures, in principle, that the results are independent of such effects. At the 
first order of the dispersion relation, we fit the evolution of the time lags as a function of z: 



<At>= aKi{z) + h{l + z), 



(5) 



where 



'1 + z) dz 



and where parameters a and b stand for extrinsic (Quantum Gravity) and intrinsic effects 
respectively. Provided t he cha nge in the defin ition of Kj. this fo rmulation is the same as the 
one used by lEUis et al.l (120061 ) . It differs from lEUis et al.l (120031 ) with respect to b parameter 
which is here expressed in the source frame of reference instead of the observer frame. 



In this paper, following lEUis et al.l (l2003l . l2006l ). we will apply the wavelet transform 
methods for noise removal and for high accuracy timings of sharp transients in 15 GRB light 
curves measured by the on-board FREGATE detector of the HETE-2 mission, assorted with 
redshift values given by the optical observations of their afterglows. The study of time lags 
between photons in various energy bands will allow us to constrain the Quantum Gravity 
scale in the linear string model as discussed above. The astrophysical effects detected in 
previous studies are also considered briefly. 

The layout of the paper is as follows: after a brief description of the HETE-2 experiment 
and gamma measurements in Section [21 we present in Section [3] the methods for de-noising 
and for the search for sharp transitions in the light curves. The results on Quantum Gravity 
scale determination are given in Section H] and possible effects other than those produced by 
the Quantum Gravity (systematic effects in the proposed analysis and astrophysical source 
effects) are presented in Section [51 Finally, the overall discussion of the results and their 
possible interpretations is presented in Section [6l 
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2. HETE-2 Experiment 

The High Energy Transient Explorer (HETE-2) mission is devoted to the study of 
GRBs using soft X-ray, medium X-ray, and gamma-ray instruments mounted on a com- 
pact spacecraft. HETE-2 was primarily developed and fabricated in- house at the MIT by a 
small scientific and engineering team, with m ajor hardware and software contributions from 



international partners in France and Japan (IDoty et al.l 120031 ). Contributions to software 



development were also made by scientific partners in the US, at the Los Alamos National 
Laboratory, the University of Chicago, and the University of California at Berkeley. Opera- 
tion of the HETE satellite and its science instruments, along with a dedi cated tracking an d 



data telemetry network, is carried out by the HETE Science Team itself (jCrew et al.ll2003l ). 
The spacecraft was successfully launched into equatorial orbit on 9 October 2000, and has 
operated in space during six years. The GRB detection and localisation system on HETE 
consists of three complementary instruments: the French Gamma Telescope (FREGATE), 
the Wide field X-ray Monitor (WXM), and the Soft X-ray Camera (SXC). T he manner in 



which the three HETE science instruments operate cooperatively is described in lRicker et al. 



(]2003f ). 



Since this study is base d on the pho t on ta gged data recorded by FREGATE, we now 



describe this instrument (see lAtteia et al.l (120031 ) for more details). The main characteristics 



of FREGATE are given in Table [TJ The instrument, which was developed by the Centre 
d'Etude Spatiale des Rayonnements (Toulouse, France), consists of four co-aligned cleaved 
Nal(Tl) scintillators, optimally sensitive in the 6 to 400 keV energy band, and one electronics 
box. Each detector has its own analog and digital electronics. The analog electronics con- 
tains a discriminator circuit with four adjustable channels and a 14-bit PHA [Pulse Height 
Analyser) whose output is regrouped into 256 evenly-spaced energy channels (approximately 
0.8 keV or 3.2 keV wide). The (dead) time needed to encode the energy of each photon is 
17 /xs for the PHA and 9 /is for the discriminator. The digital electronics processes the 
individual pulses, to generate the following data products for each detector: 

• Time histories in four energy channels, with a temporal resolution of 160 ms, 

• 128-channel energy spectra spanning the range 6-400 keV, with a temporal resolution 
of 5.24 s, 

• A "burst buffer" containing the most recent 65 536 photons tagged in time and in 
energy. 



The burst buffers are only read when a trigger occurs. The four burst buffers contain 
a total of 256k photons (64k per detector) tagged in time (with a resolution of 6.4 /is) and 
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in energy (256 energy channels). These data (also called photon tagged data) allow detailed 
studies of the spectro-temporal evolution of bright GRBs. Given the small effective area of 
each detector (40 cm^), the size of the burst buffers is usually sufficient to record all the 
ORB photons. One exception is GRB 020813 which was made of two main peaks separated 
by 60 s, and for which the burst buffers cover only the first peak. This limitation does not 
affect the analysis presented here. 

3. Description of the analysis method 

The analysis of the 15 GRBs with measured redshifts follows the steps described below: 

• determination of the time interval to be studied between start and end of burst. It is 
defined by a cut above the background measured outside of the burst region, 

• choice of the two energy bands for the time lag calculations, later called energy scenario, 
by assigning the individual photons to each energy band. The study of various scenarios 
is allowed by the use of tagged photon data provided by FREGATE for each GRB, 

• de-noising of the light curves by a Discrete Wavelet Transform and pre-selection of 
data in the studied time interval for each GRB and each energy band, 

• search for the rapid variations (spikes) in the light curves for all energy bands using a 
Continuous Wavelet Transform. The result of this step is a list of minima and maxima 
candidates, along with a coefficient characterising their regularity (Lipschitz coefficient 
a and its error (5a), 

• association in pairs of the minima and of the maxima, which fulfill the conditions 
derived from studies of the Lipschitz coefficient. 

As a result, a set of associated pairs is produced for each GRB and each energy scenario. 
The average time lag of each GRB, < At >, is then calculated and used later in the study 
of the Quantum Gravity model described in Section [H Finally, the evolution of the time 
lags as a function of Ki variable allows us to constrain the minimal value of the Quantum 
Gravity scale Eqq. 
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3.1. Use of Wavelet transforms in the analysis 

Wavelet analysi s is being increasin gly used in different fields like biology, computer 
science and physics ( iDremin et al.ll200ll ). 



Unlike Fourier transform, wavelet analysis is well adapted to the study of non-stationary 
signals, i.e. signals for which the frequency changes in time. 



These methods will be applied to our data and illustrated through figures in Section [3^ 



3.1.1. Wavelet Shrinkage 



As explained by iMallatI (119991 ). the Wavelet Shrinkage is a simple and efficient method 
to remove the noise. The Discrete Wavelet Transform (DWT) is the decomposition of a 
signal at a given resolution level L on an orthonormal wavelet basis. Such a basis can be 
defined by: 

1 ^ ft-2^n 



23 V 2J 



. = -=V' — — , (7) 



(j,n)eZ2 



where if) is the mother wavelet. The decomposition provides some large coefficients corre- 
sponding to large variations (the signal), and small coefficients due to small variations (the 
noise). Then, applying a threshold to the wavelet coefficients and performing the inverse 
transform removes the noise from the signal. In the following, the wavelet Symmlet-10 is 
used. 

There are different ways to apply a threshold to the coefficients. In this study, we first 
shift the wavelet coefficients at fine scale so that their median value is set to unity and we 
apply the so-called soft thresholding p defined by: 

X — T, ii X >T 
Prix) = { x + T, ii x<-T (8) 
0, if |x| < T. 

where x represents the wavelet coefficient and where T is a given threshold. 

The parameter T is usually chosen so that there is a high probability to be above the 



maximum level of the noise coefficients. As proposed by iDonoho fc Johnstond (119941 ). we 
use the following relation: 

T = a^/2hgN, (9) 
where a is the noise and N the number of bins of the signal. 
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It is possible to estimat e the noise value a by using the wavelet coefficients at fine scale 
( iDonoho fc Johnstond Il994j ). After ordering the N/2 wavelet coefficients at fine scale, the 
median Mx (rank is calculated. The estimator of a is then given by: 



Mx 
0.6745' 



(10) 



3.1.2. Wavelet Modulus Maxima 

The Continuous Wavelet Transform (CWT) can be used to measure the local variations 
of a signal and thus its local regularity. If a finite energy function / is considered, its CWT 
is given by: 

Wf{u, s) = f(^)^r dt, (11) 

where ip* is the wavelet function, scaled by s and shifted by u. It is possible to demonstrate 
that there cannot be a singularity of the signal without an extremum of its wavelet transform. 
A modulus maxima line is a set of points {u, s) for which the modulus of the wavelet transform 
\Wf{u,s)\ is maximum. 

It is, then, possible to detect with high precision an extremum in the data by looking 
at wavelet transform local maxima with decreasing scale, i.e. in the region of zoom of signal 
details. 

When using discrete signals, the minimum scale for looking at details should not be 
smaller than the width of one bin. In addition, the scale should not be larger than the width 
of the range in which the data is defined. So if the data is defined on the range [0, 1] by N 
bins, one must have 

A^-^ < s < 1. (12) 

Here it is important to choose a wavelet for which modulus maxima lines are continuous 
when the scale decreases. This ensures that each extremum of the light curve is localised 
by only one continuous wavelet maxima line. In the following, the second derivative of a 
gaussian, known as Mexican Hat wavelet is used. 



Like in lEllis et al.l (120031 ). extrema localised by the CWT are characterised by using the 
Lipschitz regularity, which is a measurement of the signal local fiuctuations. The Lipschitz 
regularity condition is defined as follows: / is pointwise Lipschitz a at u, if there exists a 
polynomial of maximum degree a > such that: 

\f{t)-p.{t)\<jc\t-ur, (13) 
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where /C is a constant. This mathematical expression is not used as it is in our analysis, 
instead an estimation of the Lipschitz coefficient a is obtained from a study of the decrease of 
the wav elet coefficien ts along a maxima line at fine scales. For s < Sq, the following relation 
is used jMallatlll999l ): 



iog2|iy/Ks) 



a 



- ) log2 s + k, 



(14) 



where k is a constant, sq is chosen so that it is smaller than the distance between two 
consecutive extrema. The values of the coefficient a and of its error 6a are determined by a 
fit to Equation [1] 



3.2. Data sample: light curves 

From September 2001 to May 2006, FREGATE observed 15 GRBs for which redshift 
determination was possible and for which photon tagged data are available. These GRBs 
are located from z = 0.16 to ^ = 3.37. Table [2] shows the redshift, corresponding value Ki{z) 
(see Equation [6]), duration and luminosity of each burst. 

The light curves of the 15 GRBs are shown in Figure [TJ It appears clearly that the 
signal-to-noise ratio decreases for large redshifts. The large variety of light curve shapes is a 
characteristic feature of GRBs. Different binnings are chosen depending on the duration of 
the burst. For example, we use bins of 39 ms for GRB 020813 and 12 ms for GRB 021211. 
The pre-burst sections of the light curves are used to measure the background mean and 
variance o"bck- 

After de-noising (described in Section I3.1.ip and background removal, a time range in 
which the signal is greater than abck is obtained (see top panel of Figure [2]). For GRB 030323, 
GRB 030429 and GRB 060526, a level of 0.5(Tbck is used because of a significantly lower 
signal-to-noise ratio (see bottom panel of Figure [2]). In the following, only the part of the 
light curves located in this interval are considered for measurement of the time lags. Thus, 
most of the considered extrema have a source origin. 



3.3. Choice of energy bands 

FREGATE measures the photon energies between 6 and 400 keV. After de-noising 
the light curves, extrema identification and pair association are done in energy bands. The 
selection of various energy bands is done by using the part of the spectra where the acceptance 
corrections are well understood and do not vary from burst to burst. Another important 
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consideration concerns the maximisation of the energy difference between the low and high 
energy bands. Note that the choice of contiguous energy bands provides more extrema and 
pair candidates but a smaller energy lever arm. In the following, the choice of a pair of 
energy bands will be quoted as energy scenarios. Some scenarios will produce correlated 
results due to the overlap of their energy limits. The different choices of the energy scenarios 
are summarised in Table [31 



3.4. Procedure for pair association 



Two light curves are obtained for two different energy bands, according to the method 
described in Se ction l3.3l (see Figure [3]). Noise is removed using the MatLab package WaveLab 
( lDonohdll999l ) with the wavel et shrinkage procedure described in Section [3. 1.11 (see FigureH]). 
Then the software Last Wave (lBacryll2004l ) is used to compute the Continuous Wavelet Trans- 
form (see Figure ED . It gives a list of extrema for each light curve. 

Minima and maxima are sorted in two separate groups. As maxima correspond to a 
peak of photon emission and minima correspond to a lack of photons, they do not have the 
same behavior, a priori. 



As described in Section 13.1.21 the Lipschitz coefficient a and its error 5a are deduced 
from wavelet coefficient decrease at fine scales. The derivative of the light curve f{t) at the 
position of each extremum is also computed. 

Figure E] shows that most of the extrema found by the CWT have non-zero derivatives. 
These curves can be described by a Gaussian curve centered at zero and a rather flat back- 
ground centered on large positive (negative) values for the maxima (minima) candidates. 
The width of the Gaussian curve reflects measurement errors on extrema position determi- 
nation. The fiat background corresponds to fake extrema found by our procedure. Based 
on these distributions, the following cut was applied to ensure low contribution from fake 
extrema: 

A/ 



At 



< 0.2. 



(15) 



This condition preserves all 15 GRBs in our analysis. 



A "pair" is made up of one extremum in the low energy band and one in the high energy 
band. Two values of time (ti, t2), two values of the Lipschitz coefficient (ai, ^2) and two 
values of its error {6ai, 6a2) are associated to each pair. Indices 1 and 2 are used for the 
low and high energy band, respectively. 



To build a set of extrema pairs in the most unbiased way, three variables At, Aa and 
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6{Aa) are used: 



At = t2 - h 

Aa = \a2 — Oil 



(16) 



6{Aa) = \/ba\ + ba\. 



As we focus on QG effects wliich give only small time lags, we first select possible pairs 
with |At| < 150 ms only. At this point each extrema can be used more than once depending 
on the distance between two consecutive irregularities. Then, based on the distributions of 
Aa and 5(Aq;) shown in Figures [7] and [HI the following selections are applied: 



A small value of Aa ensures that the two associated extrema are of the same kind in 
the sense of the Lipschitz regularity. The cut on 5(Aa) allows us to keep extrema mainly 
from the Gaussian peak in the 5a distribution. These cut values are valid for all energy band 
choices. 

After the selections, some extrema are used in more than one pair. To remove degeneracy 
in pair association, pairs which have the lowest At are selected. This degeneracy concerns 
only ~ 6% of the total number of pairs before the cuts (Equation [T7]) are applied. 

Table H] shows the number of pairs found for all GRBs before and after all cuts. One 
may notice the stability of the cut efficiency leading to a mean value of ~67%. 



To study the model of Quantum Gravity described in Section [T], the evolution of the 
mean time lags as a function of z was determined using the set of 15 GRBs, as illustrated 
in Figure [9] for maxima, minima and combination of both for scenario #2. We show on the 
same figure the fit of the data points with Equation [51 The results (a, 6) of the fits for all 
scenarios are summarised in Table [51 As explained in Section [H the parameter a depends 
on the Quantum Gravity scale while the parameter h refiects intrinsic source effects. Both 
parameters were found to be strongly correlated, as shown in Figure [TUl which represents 
the 95% Confidence Level (CL) contours for a and h. In spite of the fact that maxima and 
minima present similar exclusion domains, it should be noticed that most ellipse centres are 
grouped around zero value for a and h in case of maxima and when both minima and maxima 
are considered, whereas they are slightly shifted towards values of a > and 6 < for the 




(17) 



4. Results on Quantum Gravity scale 
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minima. However, the fit results suggest no variation above ±3 a, so that in the following, 
we derive the 95% CL lower Quantum Gravity scale limit, assuming no signal is observed. 

The dependence of the Quantum Gravity scale parameter Eqg on z may also be con- 
st rained__by_a direct study of the sensitivity of the 15 GRB data to the model proposed 



by lEllis et al.l (120001 ). We have built a likelihood function following the formula: 



2 

where M is the energy and the x^(M) is expressed as 



e.p(-^l, (18) 



i=l ^ 

where the index i corresponds to each GRB and Nqrb < 15 is the total number of GRBs 
with at least one pair, for a given scenario. The dependence of a on M as predicted by the 
considered model of Quantum Gravity is given by: 

The mean values of energy are computed for each GRB and each scenario from the relation: 

A<^> = <E>2-<^>i, (21) 

where indices 1 and 2 represent the low and the high energy band, respectively. The averaged 
values of A < E' > for all bursts are given in Table [3] for each scenario. In this study, a 
universality of the intrinsic source time lags has been assumed. The average value h (and its 
error o"^) is obtained as the weighted mean of the values hk (and their errors ak) from the 
previous two-parameter linear fit: 

h = ^ and a-, = , (22) 

where the index k corresponds to each scenario and Wk = l/o"^- Using the values of bk and 
ak given in Table E, one gets b = 0.0023 ± 0.0026 for maxima, b = -0.0282 ± 0.0038 for 
minima and b = —0.0069 ± 0.0035 for both minima and maxima. 

Figure [TT] presents the evolution of x^(^)/ndf around its minimum x^j„(M)/ndf for 
maxima, minima and all extrema. All scenarios fulfill the condition x^^„/ndf < 2. Like in 
the two-parameter linear fit, these curves show a different behaviour in case of the maxima 
and minima: no significant preference of any value of M is observed for most of the scenarios 
for the maxima, whereas a preferred minimum seems to be found for the minima. 
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The 95% CL lower limit on the Quantum Gravity scale is set by requiring A(x^/ndf) 
to vary by 3.84 from the minimum of the function. The values of limits on Eqg (shown 
in Table E]) are obtained for the 14 energy scenarios. In good agreement with the slope 
parameter a, all values for minima and maxima are within the lO^'^-lO^^ GeV range. No 
correlation with energy band choice or with any other cut parameter is found. As there is 
no reason to choose any particular value of the lower limit on the Quantum Gravity scale 
obtained for the 14 scenarios, we choose scenario #3, which provides the strongest constraint 
on the QG scale. This gives us a lower limit on Eqg of 3.2 x 10^^ GeV for the maxima and 
a smaller value of 7.5 x 10^^ GeV for the minima. Combining maxima and minima, the best 
lower limit for the QG scale remains of 2 x 10^^ GeV. 

5. Discussion on various origins of time lags 

5.1. Measurement uncertainties in the time lag determination 

The time lag determination procedure proposed in this analysis is subject to various 
systematic effects due to: 

• the experimental precision on time measurement and energy estimation by the FRE- 
GATE detector, 

• the biases that could be introduced by the use of the wavelet methods, 

• the procedures employed for the extrema selection and pair association. 

The precision on the individual photon time measurement is governed by the precision 
on the HETE-2 internal clock time synchronised to UTC by a spatial GPS whithin few /xs, 
resulting in an overall precision value below 1 ms, well below the bin width in the CWT 
decomposition. The energy resolution of the FREGATE detector in the studied domain does 
not exceed 20% and has little impact on the results of the present analysis. The error on 
redshift determination by the optical observations of the GRBs, can be considered negligible 
as well. 

The systematic effects of the wavelet decomposition may have impact not only on the 
de-noising of the light curves but also on the extrema localisation, producing different biases 
in the maxima and minima identification. This could be the explanation of the different 
behavior of curves for maxima and minima. Since systematic effects may be different for 
maxima and minima, we consider the result obtained by the combination of both maxima 
and minima as more reliable. 
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The result of the de-noising procedure depends on two parameters: the type of the 
wavelet function and the level of the decomposition L. Wavelet functions Symmlet-10 and 
Daubechies-10 were considered. Both have ten vanishing moments but Symmlet is more 
symmetric. Both wavelets give comparable r esults: no discrep ancy was observed between 



extrema positions. Following the approach by lEUis et al.l (120031 ). the Symmlet-10 wavelet is 



used for the results given in this paper. As far as the decomposition level L is concerned, 
different values were tested. The lower the value of L is, the smoother is the obtained light 
curve and the fewer extrema are found. The value L = 6 was chosen because it allows us to 
keep a significant level of detail without adding too much noise. The noise fluctuation cut at 
la above background was also raised to higher values, consequently decreasing the number 
of extrema candidates. However, no significant change in the results was observed. The cuts 
applied on values of the derivative at each extrema found with the CWT were studied. A 
more stringent cut of < 0.1 value was also applied to reject the fake extrema. This cut rejects 
about 40% of pairs giving less significant results, while the cut value of 0.2 (Equation [15]) 
rejects about 15% of pairs. Concerning the cuts on Aa and 6{Aa), two selections, more 
severe than Equation [TTl have been investigated: 

Aa < 0.2 j Aa< 0.4 

(5(Aa) < 0.045 \ 5(A«) < 0.02. ^ ^ 

The different choices of cuts on the extrema selections and pair association had little 
effect on the final limits, even if the statistical sensitivity of the studied sample was reduced. 
In particular, fewer extrema candidates were found in the light curves, leading to a smaller 
GRB sample. 



5.2. Intrinsic source effects 



Even in the restricted range of energies in our study, the GRB light-curves are not 
"perfect" signals, since they exhibit intrinsic time lags between high and low energies, which 
vary from burst to burst. It has been known for a long time, that the peaks of the emission 



are shorter and arrive earlier at higher energies (IFenimore et al.l 1 19951 : iNorris et al.l Il996 



and references therein). These intrinsic lags, which have a sign opposite to the sign expected 
from Lorentz violation, have a broad dispersion of durations, complicating the detection of 
Lorentz violation effects. Therefore the effect of Lorentz violation must be searched with a 
statistical study analysing the average dependence of the lags with z, and not with a single 
GRB. 



In particular, a strong anticorrelation of spectral lags with luminosity has been found 
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by iNorris et al.l (120001 ) . At low redshifts, we detect bright and faint GRBs which present a 
broad distribution of instrinsic lags whereas at high redshifts, only bright GRBs with small 
intrinsic lags are detected. This effect could mimick the effect of Lorentz violation due to 
the non uniform distribution in luminosity of the GRBs in our sample. 



Following lEUis et al.l (120061 ). we use the generic term of "source effects" for the intrinsic 
lags discussed above. Source effects must be carefully taken into account if one wants to 
derive meaningful limits on the magnitude of Lorentz violation. Indeed , the source effect s 
mentioned afore could explain the positive signal recently reported by lEUis et al.l (120061 ). 
There are at least two ways to take into account the source effects: 



• Their modelisation, which allows to understand the impact of the sample selection on 
the final result. 

• The selection of a GRB sample homogeneous in luminosity, which minimises the impact 
of source effects. 



The modelisation of the source effects is beyond the scope of this paper (but it will 
become more and more important in future studies, when the increasing number of GRBs 
will allow to place stronger constraints on Quantum Gravity effects). Here, with limited 
statistics, we made a test by performing our analysis on a restricted GRB sample, almost 
homogeneous in luminosity. The HETE-2 sample we use comprises 15 GRBs with redshifts 
ranging between 0.1 and 3.4, biased with respect to the luminosity population. The high 
redshift part of the sample is mainly populated by the high luminosity GRBs for which the 
source time lag effects are expected to be the weakest. The independent analysis of a GRB 
sample with luminosity values L51 > 8 erg s~^ (10 GRBs) provided similar results to those 
obtained with full sample of 15 GRBs. The sensitivity was lower because of a decreased 
statistical power of the restricted sample and a smaller lever arm in the redshift values. 



6. Summary 

We have presented in this paper the analysis of the time structure of 15 GRB light curves 
collected by the HET15-2 mission in years 2001-2006, for which the redshift values have been 
measured. This sample was used for time lag measurements as a function of redshift. These 
time lags can originate either from the astrophysical sources themselves or from a possible 
Quantum Gravity signature. The latter case would result in energy and redshift dependence 
of the arrival time of photons. 
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We have used one of the most precise methods based on the wavelet transforms for the 
de-noising of the hght curves and for the locahsation of sharp transitions in various energy 
ranges. The maxima and the minima in the hght curves have been studied separately as 
well as together. In particular, the dependence on the Quantum Gravity scale for the 
maxima shows a continuous decrease with energy scale parameter M for most of the energy 
scenarios, whereas a minimum value may be detected for the minima. For 14 choices of 
the energy bands, the observed slope as a function of redshift for the maxima does not 
exceed 3 a variation from a zero value. The preferred value in case of the minima is situated 
between lO^'' and 10^^ GeV for almost all considered energy scenarios. We use minima and 
maxima together to determine the 95% CL lower limit on the Quantum Gravity scale for 
two reasons. The first reason is that we cannot exclude zero values at 95% CL for a and 
b parameters of the two parameter fit reflecting QG and source effects. The second is that 
there is only a slight difference between the results for minima and maxima. In most of 
the studied energy scenarios, the lower limit value is of the order of 10^^ GeV and can be 
considered as competitive considering the modest energy gap of ~130 keV provided by the 
IIETE-2 data. In this study, we do not correct for the astrophysical source effects, except in 
case of the analysis performed on a restricted GRB sample, homogeneous in luminosity. 

The impact of the selections and cuts on the obtained results has also been investigated. 
The most important contribution to the systematic effects on the background suppression 
comes from the decomposition level choice in the discrete wavelet transformation used in 
the de-noising procedure. All cuts applied in the extrema identification and in the pair 
association have also been varied, and lead to compatible results with those obtained from 
the optimised selections. The other important factor in this analysis is related to the energy 
lever arm in the selected energy scenarios. Here, we would like to underline the importance 
of the access to the full information on energy and time variables provided by the FREGATE 
detector, delivered on an individual photon basis. The FREGATE precision on photon time 
measurement is also a crucial parameter for the proposed analysis. 

In summary, studies of time lags from the position of all the extrema in the light curves 
of the HETE-2 GRBs allow us to set a lower limit on the Quantum Gravity scale of 

Eqg > 2 X 10^^ GeV. 

The curves for the minima of the light curves show a prefered minimum between 10^^ and 
10^^ GeV. All studies of the systematic effects do not change these results in a significant 
way. Further improvement of the limits on the Quantum Gravity energy scale would need 
an order of magnitude larger statistics and a larger energy lever arm. The next step in this 
direction, would be th e analysis of an enlarg ed sample of HETE-2 GRBs, with measured 



pseudo-redshift values (jPelangeon et al.ll2006l ) 
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A. Time lags and cosmological effects 

The relation between time and redshift is given by 

dt = -Ho^ (Al) 
(1 + z) h[z) 

where h{z) is defined by Equation H] and where Hq is the Hubble constant. 
During a time dt, a particle with velocity u travels a distance of: 

XL dz 

dl = udt = -Hn ^ 7 rvT^- (A2) 

" {l + z)h{z) 

It is important to note here that this distance is measured at redshift z. The same distance 
measured at the redshift of the oberver is then given by 

dh = -Ho ^ (A3) 
So, two particles with velocities different by Am travel different distances by: 

From Equation [2] (Section [1]), we deduce that two photons with an energy difference AE 
at redshift present a velocity difference at redshift z of: 

A„=-c^^il±il. (A5) 
Eqg 
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The integration of equation IA4I provides the final formula (Equation [3] in the text) for 
the time lag when the two photons are emitted at the same time: 
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Table 1. Main parameters of the FREGATE detector. 



Parameter Value 

Energy range 6—400 keV 

Effective area (4 detectors, on axis) 160 cm^ 

Field of view (HWZM) 70° 

Sensitivity (50-300 keV) lO"''' erg cm-^ s'^ 

Dead time 10 fis 

Time resolution 6.4 /is 

Maximum acceptable photon flux 10* ph cm~^ 

Spectral resolution at 662 keV ~8% 

Spectral resolution at 122 keV ~12% 

Spectral resolution at 6 keV ~42% 



Table 2. 15 GRBs observed by HETE-2 from September 2001 to May 2006. 



GRB z Ki T90 (s) L6i(ergs-i) 



GRB 


050709 


0, 


.16 


0, 


.17 


0.1 




GRB 


020819 


0, 


.41 


0, 


.44 


12.6 


0.64 


GRB 


010921 


0, 


.45 


0. 


.49 


21.1 


1.31 


GRB 


041006 


0, 


.71 


0. 


.79 


19.0 


5.46 


GRB 


030528 





.78 


0, 


.87 


21.6 


1.21 


GRB 


040924 


0, 


.86 


0, 


.96 


2.7 


9.10 


GRB 


021211 


1, 


.01 


1, 


.13 


2.4 


11.97 


GRB 


050408 


1, 


.24 


1, 


.38 


15.3 


9.51 


GRB 


020813 


1, 


.25 


1. 


.40 


89.3 


33.51 


GRB 


060124 


2, 


.30 


2, 


.49 


18.6 


43.44 


GRB 


021004 


2, 


.32 


2, 


.51 


53.2 


9.28 


GRB 


030429 


2, 


.65 


2, 


.83 


10.3 


11.24 


GRB 


020124 


3, 


.20 


3, 


.32 


46.4 


53.52 


GRB 


060526 


3, 


.22 


3, 


.34 


6.7 


21.21 


GRB 


030323 


3, 


.37 


3. 


.47 


27.8 


11.92 



Note. — T90 is defined as the time during which 5% 
to 95% of the total observed counts have been detected in 
the energy range 30—400 keV. The luminosity is expressed in 
units of 10^1 ergs-l . Luminosity of GRB 050709 is not given 
because this burst is much shorter than the other bursts. 
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Table 3. Different choices of the energy bands studied in HETE-2 GRB analysis, quoted 
as energy scenarios in the text. Values of < AE >, averages for all GRBs are given in the 

third column. 



Scenario 


Energy band 1 


Energy band 2 


Mean < AE > 


#1 


20-35 keV 


60-350 keV 


117.6 kcV 


#2 


8-30 keV 


60-350 keV 


127.2 kcV 


#3 


8-20 keV 


60-350 keV 


130.2 keV 


#4 


8-20 keV 


30-350 keV 


85.0 keV 


#5 


8-30 keV 


30-350 keV 


82.0 keV 


#6 


8-20 keV 


40-350 keV 


102.8 keV 


#7 


8 30 kcV 


40-350 keV 


99.8 keV 


#8 


8-40 kcV 


40-350 keV 


97.9 keV 


#9 


20-35 keV 


40-350 keV 


90.1 keV 


#10 


8-20 keV 


50-350 keV 


116.9 keV 


#11 


8-30 keV 


50-350 keV 


113.9 keV 


#12 


8-40 keV 


50-350 keV 


112.0 kcV 


#13 


8-50 keV 


50-350 keV 


110.4 keV 


#14 


20-35 keV 


50-350 keV 


104.2 keV 



Table 4. Number of pairs found for all GRBs and each energy scenario. 



Scenario 


Before cuts 


After cuts 


Efficiency 


#1 


165 


118 


72% 


#2 


118 


88 


75% 


#3 


139 


97 


70% 


#4 


132 


93 


70% 


#5 


109 


77 


71% 


#6 


141 


91 


65% 


#7 


114 


73 


64% 


#8 


109 


65 


60% 


#9 


159 


102 


64% 


#10 


145 


90 


62% 


#11 


122 


79 


65% 


#12 


112 


72 


64% 


#13 


103 


61 


59% 


#14 


159 


111 


70% 




Mean efficiency : 


67% 



Note. — Number of pairs is obtained for both maxima 
and minima. 
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Table 5. Results of the fit < At >= a Ki{z) + h {1 + z) for all scenarios. 





Minima 


Maxima 


All 


Scenario 


a 


b 


a 


6 


a 


b 


#1 


0.0027±0.0116 


0.0001±0.0070 


-0.0780±0.0426 


0.0613±0.0279 


-0.0422±0.0347 


0.03r59±0.0224 


#2 


0.0487±0.0288 


-0.0317±0.0184 


0.0525±0.0503 


-0.0497±0.0327 


0.0400±0.0373 


-0.0371±0.0240 


#3 


0.0481±0.0191 


-0.0318±0.0116 


-0.0034±0.0095 


0.0048±0.0065 


0.0169±0.0086 


-0.0089±0.0059 


#4 


0.1144±0.0433 


-0.0917±0.0272 


0.0807±0.0378 


-0.0645±0.0246 


0.0631±0.0324 


-0.0493±0.0218 


#5 


0.0447±0.0296 


-0.0245±0.0182 


0.0274±0.0530 


-0.0330±0.0312 


0.0409±0.0477 


-0.0345±0.0299 


#6 


0.0812±0.0130 


-0.0527±0.0076 


-0.0102±0.0081 


0.0101±0.0051 


0.0211±0.0142 


-0.0105±0.0101 


#7 


0.0164±0.0527 


0.0146±0.0308 


0.0225±0.0169 


-0.0173±0.0117 


0.0009±0.0383 


0.0105±0.0243 


#8 


0.1056±0.0567 


-0.0449±0.0308 


-0.0283±0.0160 


0.0122±0.0095 


0.0159±0.0142 


-0.0139±0.0100 


#9 


0.0651±0.0327 


-0.0524it0.0187 


0.0286±0.0406 


-0.0272it0.0246 


-0.0044it0.0149 


0.0058±0.0097 


#10 


0.0756±0.0330 


-0.0338±0.0188 


-0.0011±0.0089 


0.0080±0.0059 


0.0139±0.0296 


0.0009±0.0195 


#11 


0.0493±0.0274 


-0.0263±0.0179 


0.0341±0.0251 


-0.0242±0.0163 


0.0584±0.0322 


-0.0389±0.0220 


#12 


0.0353±0.0305 


-0.0246±0.0198 


-0.0373±0.0308 


0.0379±0.0183 


0.0209±0.0207 


-0.0136±0.0142 


#13 


0.1001±0.0493 


-0.0417±0.0271 


0.0313±0.0126 


-0.0206±0.0078 


-0.0366±0.0432 


0.0348±0.0280 


#14 


0.0604±0.0525 


-0.0117±0.0301 


-0.0444±0.0385 


0.0347±0.0257 


-0.0358±0.0286 


0.0310±0.0188 



Note. — Errors are obtained from the fit and normalised to jndj ~ 1. 



Table 6. Obtained 95% CL lower limit on Eqg (GeV). 



Scenario 


Minima 


Maxima 


All 


#1 


5.9 X 


10" 


5.8 X 10" 


5.7 X 10" 


#2 


5.8 X 


10" 


7.8 X 10" 


9.0 X 10" 


#3 


7.5 X 


10" 


3.2 X 10^5 


2.0 X 10^5 


#4 


3.6 X 


10" 


5.8 X 10" 


6.3 X 10" 


#5 


3.8 X 


10" 


4.6 X 10" 


4.6 X 10" 


#6 


4.7 X 


10" 


1.8 X 10^5 


1.3 X IQis 


#7 


2.3 X 


10" 


1.7 X 10^5 


4.5 X 10" 


#8 


2.1 X 


10" 


1.2 X IQis 


1.8 X IQis 


#9 


4.0 X 


10" 


5.8 X 10" 


1.0 X 10^5 


#10 


3.6 X 


10" 


1.9 X 10^5 


7.3 X 10" 


#11 


5.3 X 


10" 


8.6 X 10" 


7.9 X 10" 


#12 


5.8 X 


10" 


5.0 X 10" 


1.2 X 10^5 


#13 


2.4 X 


10" 


1.2 X IQis 


5.2 X 10" 


#14 


2.3 X 


10" 


7.7 X 10" 


6.8 X 10" 



-24- 




Fig. 1. — Light curves of the 15 GRBs detected by HETE-2 in the energy range 6-400 keV. 
The bursts are sorted by increasing z from top left to bottom right. X-axes are graduated 
in seconds. 
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Fig. 2. — Noisc-frcc, background subtracted light curve of GRB 041006 (top) and of 
GRB 030323 (bottom). For GRB 041006, the horizontal line corresponds to Idbck- For 
GRB 030323, the horizontal line corresponds to 0.5cTbck- The hatched areas show the parts 
of the light curve that do not contribute to the analysis. 
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Fig. 3.— Light curves of GRB 041006 for the energy bands 8-30 keV and 60-350 keV. 
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Fig. 4. — De-noised light curves of GRB 041006 for the energy bands 8-30 keV and 
60-350 keV. These curves were obtained with the wavelet Symmlet-10 with the soft thresh- 
olding proceedure at the level L — 6. 
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Fig. 5. — The Continuous Wavelet Transform of the hght curve between 8 and 30 keV of 
GRB 041006. The top panel shows the denoised data. The middle panel shows the wavelet 
transform, i.e. the value of \Wf{u, s) \ for each [u, s). The horizontal and vertical axes give 
respectively u and a = log2 s. The bottom panel shows the modulus maxima of Wf{u, s). 
Each maxima line corresponds to a singularity in the signal. 
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Fig. 6. — Distribution of the derivatives of the hght curve for maxima (top) and minima 
(bottom) located at to. Only extrema in the range found in Section 13.21 are considered. 
The data is fitted for —0.4 < {Af /At)t=to — ^-^ with the sum of two Gaussian curves. The 
contribution of fake extrema is negligible for |A//At| < 0.2. 
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Fig. 7. — Distribution of Aa for all candidate pairs in scenario #2. Cut is shown by an 
arrow. 
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Fig. 8. — Distribution of (5(Aq!) in the scenario #2. The cut is based on the position of the 
peak, determined by a Gaussian fit. We set the cut at three times the mean value of (5(Aq;). 
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Fig. 9. — At in seconds, as a function of z for scenario #2, considering maxima only (top), 
minima only (middle) and both minima and maxima (bottom). The curves show the result 
of a fit to Equation [51 
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Fig. 10. — 95% CL contours for a and b from the two-parameter fit for the fourteen scenarios, 
for maxima only (top), minima only (middle) and both minima and maxima (bottom). The 
boxes at the bottom left of the plots show the position of contour centers. 
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Fig. 11. — Evolution of function of M, for maxima only (top), minima only (middle) and 
both minima and maxima (bottom). 



